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Abstract 

The white light during M dwarf flares has long been known to exhibit the 
broadband shape of a T w 10 000 K blackbody, and the white light in solar 
flares is thought to arise primarily from Hydrogen recombination. Yet, a current 
lack of broad wavelength coverage solar-flare spectra in the optical/near-UV 
prohibits a direct comparison of the continuum properties to determine if they 
are indeed so different. New spectroscopic observations of a secondary flare 
during the decay of a megaflare on the dM4.5e star YZ CMi have revealed 
multiple components in the white-light continuum of stellar flares, including 
both a blackbody-likc spectrum and a hydrogen recombination spectrum. One 
of the most surprising findings is that these two components are anti-corrclated 
in their temporal evolution. We combine initial phenomenological modeling of 
the continuum components with spectra from radiative-hydrodynamic models 
to show that continuum veiling gives rise to the measured anti-correlation. This 
modeling allows us to use the components' inferred properties to predict how 
a similar spatially resolved, multiple-component white-light continuum might 
appear using analogies to several solar flare phenomena. We also compare the 
properties of the optical stellar flare white light to Ellerman bombs on the Sun. 

Keywords: white light flares, solar-stellar connection, radiative transfer, Eller- 
man bombs 



1. Introduction 

In both solar and stellar flares, the near-UV and optical (white-light) continuum 
emission is an energetically important but unexplained phenomenon. On the Sun, 
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the white- light continuum appears in small regions of transien t emission that are 
spatially and tempo r ally coincident with hard X-ray bursts (IRust and Hegwei , 
19751 : iHudson et all 1 19921 : iNeidig and Kanel Il993t iFletcher et all 120071) . This 



relation suggests that the origin of the white light is related to the energy 
deposited in the lower atmosphere by nonthermal electrons accelerated during 
flares. Broad wavelength coverage spectral observatio ns are sparse and date back 
to several large sol a r flares from the 1970s and 1980s (IMachado and Rusu Il974 ; 
ffiel Il982t iNeidid . Il983l iDonati-Falchi. Smaldone. and FalcianiL Il984h These 



spectra are consistent with continua arising primarily from the hydrogen B aimer 
continuum and H _ emission. 

Whereas the largest solar flares emit < 10 32 erg s in the white-light continuum 
and l ast not much longer than ten minutes (e.g., Neidig. Grosser, and Hrovatl . 
19941 ) . the white- light emis sion on active lower-mass M dwarfs can reach > 10 34 



ergs and persist for hours (jHawlev and Pettersenl . Il991t Kowalski et all 12010a , 
hereafter KlOa). More is known about the spectral shape of the white light 
during M dwarf flares (many spectrographs can easily obtain low-resolution, 
broad wavelength spectra of stellar flares, but it is difficult to place a spec- 
trograph slit over solar white-light kernels, which are intermitte nt and largely 
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Zhilvaev et al. . 2007|) and optical/NUV spectr a ( Hawlev and Pettersen . 



Eason et al . 19921 Garcia- Alvarez et al . 2002 : Fuhrmeister et 0^2008 



KlOa). In contrast to solar observations, the spectral shape during M dwarf 
flares suggests a hot blackbody with temperatures between 8500 — 11000 
K. Areal coverages of this component are typically < 0.1% of the visible stellar 
hemisphere, which implies a compact geometry like those observed in white light 
at the footpoints of flare arcades on the Sun. Although this blackbody component 
seems to be nearly ubiquitous during (large) stellar flares, it is not p redicted even 



by th e most recent ID radiative hydrody namic (RHP) flare models (I Alfred et al. 
20061) produced with the RADYN code (|Carlsson and Steinl . Il994 Il995l Il997t) 



The "megaflare" of UT 16 January 2009 on the dM4.5e star YZ CMi is one of 
the longest lasting and most energetic flares observed on a low-mass single star. 
Low- resolution spectra (3350 - 9260A) were obtained in the flare's decay phase, 
which was elevated between 15 and 37 times the quiescent level and contained 
many secondary peaks. More than 160 spectra were obtained over 1.3 hours, 
and simultaneous U-band photometric observations of the entire flare event were 
provided by the NMSU lm telescope. A detailed description of the observations 
and data reduction is given in KlOa. 

In KlOa, two continuum components were necessary to fit the blue (3350- 
5500A) spectra: a hydro gen Balmer continu um (BaC) component as predicted 



by the RHD models of lAllred et al 



(|2006l ) and a T « 10 000 K blackbody 
component. An intriguing anti-correlation was found between the temporal evo- 
lution of these two components: the blackbody emission increased when the BaC 
decreased, and vice versa (see Figure Id of KlOa). In this article, we revisit this 
anti-correlation and provide an explanation for it using the phenom enological 
models of the secondary flare spectra from iKowalski efall (|2010bl hereafter 
KlOb). Finally, we show how each component of this flare might appear in the 
context of a solar flare "arcade" . 
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2. Anti-Correlated Continuum Components and Continuum Veiling 

The anti-correlation between the blackbody and BaC components (KlOa) can be 
understood qualitatively using Figure [TJ where we show the spectral evolution 
of the total "flare-only" flux (denoted here as F^) during the rise phase of the 
secondary flare at t ~ 130 minutes. The spectra are color-coded to the nearest 
(in time) U-band measurement in the inset panel. At times prior to and near 
the beginning of the secondary flare (black, purple, and dark blue spectra), two 
distinct continuum components are clearly present in the spectra. The best-fit 
blackbody (short-dashed line) accounts for most continuum emission at A > 
4000A, whereas the BaC emission above the blackbody is conspicuous at A < 
3750A. During the rise and at the peak of the secondary flare (green, yellow, 
and red spectra), the BaC component seemingly disappears and the best-fit 
blackbody (long-dashed line) can fit the continuum shape throughout the entire 
wavelength range. KlOa showed that the H7 line flux exhibits an anti-correlated 
relation with the blackbody component. This effect is also present in Figure [TJ In 
the red (secondary flare peak) spectrum, the continuum at A « 4200A is highest, 
yet the peaks of the hydrogen Balmer lines are lowest. 

The secondary flares at t « 95 minutes and t « 130 minutes are events during 
which the blackbody flux becomes stronger while the BaC flux becomes weaker. 
KlOa quantified this as an increasing filling factor (areal coverage; percent of 
stellar disk) of the blackbody (with constant temperature, T = 10 000 K) during 
the rise phase of each secondary flare. In Figure [TJ we present an alternative 
interpretation. The blackbody curves (dashed lines) were fit to the spectra by 
allowing both the temperature and filling factor to vary. The best-fit blackbody 
temperatures and filling factors are T m 10 400 K and X B b ~ 0.1% (at t = 122.9 
minutes; short-dashed line) and T w 13 000 K and Xbb ~ 0.1% (at t = 130 
minutes; long-dashed line). Strikingly, if both parameters arc allowed to vary 
when fitting a blackbody function to these total flare spectra, the temperature 
increases by « 2500 K while the filling factor remains approximately constant. 
Fitting a blackbody to the total flare spectrum (either by holding T constant, 
or by allowing X and T to vary) gives only the average properties of the entire 
flaring region at that time. We next show that these interpretations can be 
improved by isolating the newly-formed flare emission. 

KlOb found that the isolated flare emission (denoted here as F£) during 
the secondary flare's rise phase resembles the spectrum of a hot star, with the 
defining features being the hydrogen Balmer continuum and lines in absorption 
and a steeply rising continuum towards the blue at A > 4000A (Figure lb of KlOb 
shows that the new flare emission, obtained by subtracting the pre-secondary 
flare spectrum (average of three black and purple spectra around t = 123.4 
minutes) from the average of two green-colored spectra around t = 126.5 minutes 
in Figure [TJ of this article, is very similar to the spectrum of the AO star Vega) . 
The observed anti-correlation between the continuum components in Figure [JJ is 
a result of a 'hot star spectrum' forming during the secondary flare. The hot star 
( "blackbody-like" ) spectrum causes an increase in the continuum at A ~ 4200A 
by an amount, -F4200J whereas an increase in the continuum at A ~ 3500 A occurs 
by only w 0.6 x F± 200 . In other words, the flux in the continuum on both sides of 
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the Balmer jump increases, but the continuum at A rj 4200A increases by a larger 
amount. Thus, the apparent decrease in the total amount of BaC in emission 
from t = 123 minutes to 130 minutes occurs as a result of 'continuum veiling' 
(similar to the c ontin u um veiling observed for accreting T Tauri stars - see, e.g., 
Hartigan et al. ( 1989h : Hessman and Guenther (1997); Herczeg and Hillenbrand! 



(2008)). 
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Figure 1. A scries of 16 flare spectra obtained prior to the onset and through the peak of 
a secondary flare. The quiescent spectrum from 24 November 2008 has been subtracted, as 
in KlOa. The U-band light curve (inset) is color-coded to the spectrum obtained closest in 
time. The best-fit blackbody curves to the black- and red-colored spectra are shown as the 
short-dashed (T Ri 10 400 K) and long-dashed lines (T ps 13 000 K), respectively. 



Combining Continuum Components using Phenomenological Hot 
Spot Models 



In KlOb, the F^ emission was modelled phen omenologically with the static 
radiative transfer code, RH (jUitenbroekl . 120011) . as a temperature bump ("hot 
spot") with peak temperature T = 20 000 K, placed near the photosphere (below 
the temperature minimum) of the quiescent M dwarf atmosphere. Here, we use a 
su m of individual hot spots and the RHD model spectrum (hereafter RHDF11) 
of Alfred et al. l|2006l) to model the total flare emission (F' x ) at two times during 
the megaflare on YZ CMi. Figure [5] shows flare spectra from Figure Q] averaged 
around t = 123 .4 minutes (F0 1 ; grey) and at t = 126.5 minutes (-FT'; black). 
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These are the spectra corresponding to times immediately before and nearly 
half- way up the rise phase of the secondary flare, respectively (i.e., the same two 
spectra presented in Figure la of KlOb but with the quiescent level subtracted). 



Figure 2. Flare spectra at i = 123.4 minutes and 126.5 minutes from Figure ^ 
arc shown in grey and black, respectively. The composite model spectra are shown 
in blue, short dashes (F^rhdfii x Xrhdfii + Fj^hsi x -Xhsi) and red, long dashes 
(^a.rhdfii x ^rhdfii +^A,HSi x ^HSl + ^A,HS2 X X H S2)- The RHDF11 spectrum is shown 
as the thin black line. The continuum veiling effect is apparent from the different heights of 
the Balmcr discontinuity at A = 3646A. 

To model FQ', we add the RHDF11 spectrum and a hot spot (HS1) with 
Imax = 12 000 K (keeping the other parameters the same as described in KlOb) 
with a ratio of filling factors of 10 : 1 and Arhdfii = 1.2%, as found in KlOa. 
The total model spectrum is shown as the light blue (short dashes) curve in 
Figure [21 As in KlOb, wc model the secondary flare as a hot spot (HS2) with 
T max = 20 000 K. Adding HS2 to RHDF11 and HS1 gives the red (long dashes) 
spectrum in Figure [2] The areal coverage of HS2 is 0.46 as large as the areal 
coverage of HS1. These model spectra match the observed continuum levels at 
all wavelengths in the figure. Moreover, the continuum veiling / anti-correlation 
effect is readily apparent: the height of the Balmer jump at A = 3646A decreases 
from black (no hot spots) to blue (one hot spot) to red (two hot spots). The 
model fluxes are slightly lower than the observed fluxes at A < 3750A. At these 
wavelengths, there is a forest of metallic lines (e.g., Fe I, Fe II) that are blended in 
our low-resolution spectra; our model is satisfactory in matching the underlying 
level which is likely closer to the actual level of the BaC. Note that in KlOa, 
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we intentionally used only the BaC from the RHDF11 spectrum to model the 
continuum, whereas the RHDF11 predictions for the Paschen continuum and 
photospheric-backwarming components are included in the total fluxes in this 
work. 

The origin of the hot spots is unknown, as they are not predicted by self- 
consistent radiative hydro dynamic models th at employ a solar-type non-thermal 
electron heating function ( Allred et al. . 2006;). In the proposed continuum model, 
we use the fewest number of components necessary to fit the overall shape and 
reproduce the anti-correlation. However, HS1 may represent a sum of individual 
decaying hot spots from previously heated flare regions (see below). We are 
working to produce a grid of phcnomcnological models which will be used to 
constrain the column mass of the hot spots, the detailed temperature evolution, 
and the uniqueness of the continuum fit. Radiative hydrodynamic models (with 
RADYN) of the gradual phase are also forthcoming and will provide a more 
accurate representation of the BaC and photosphcric backwarming in the decay 
phase. 



4. The Solar Analogy 

Figure [3] shows how a spatially resolved observation of the YZ CMi megaflare 
might have appeared. We use the continuum components and filling factors to 
make analogies to several flare structures and phenomena observed in large solar 
flare arcades. The complex morphology of the U-band light curve leads us to 
speculate that the YZ CMi megaflare involved a large arcade, or several large 
arcades of flare loops. The main features of Figure [3] are the following: 

• BaC (yellow ribbons): Throughout the spectral observations (72 minutes 
< t < 149 minutes), the hydrogen BaC and lines were highly elevated 
and decreasing (likely from the initial flare peaks in the U-band light 
curve), implying that this emission had originated from the footpoints of 
a previously heated magnetic arcade in the flaring chromosphere. These 
may manifest as a complex of flare ribbons, as is commonly observed in 



Ha during solar flares (e.g., Rust and Hegwer . 1975 : Berlicki et all 2004 : 
Balasubramaniam et all . 201Clh . In some solar flares (e.g., Neidig, 1983 ). the 



BaC appears to have a spatial morphology that is more compact than an 
extended Ha ribbon. Spectroscopic observations arc needed to compare the 
plasma properties and conditions of BaC and Balmer line emitting ribbons 
and kernels. 

HS1 (purple spots): Immediately prior to the secondary flare beginning 
at t « 123 minutes, a series of phenomenological hot spots (HS1) are 
present near the photosphere. These hot spots were formed during the 
previous secondary flares (i.e., at t m 65 minutes, ps 95 minutes; see KlOa), 
and they are emitting from a total source size that is w 1/10 as large as 
the area of the chromosphcric flare region (e.g., Ha ribbons). The spectra 
of these hot spots have the hydrogen BaC and lines in absorption. They 
might be similar to the compact white-light kernels during solar flares, as 
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i nlWang. Fang, and Ming-DeDing] (|2007l) . lFletcher et all (|2007h . I Jess et al 
( 20081 ). or thev may be similar to Ellerman-bomb phenomena (see Section 
14. ip . Also emitting from the photosphere is a larger region heated from 
chromospheric (BaC) backwarming; we assume the size of this backwarmed 
region is the same size as the flaring chromosphere. 

HS2 (white spot): The secondary flare at t « 130 minutes is the result 
of the formation of a new hot spot (HS2), hotter and smaller than HS1 but 
at the same column mass. At this time, we see a sudden decline in the BaC 
flux. When all of these components are unresolved, as in our stellar spectra, 
continuum veiling gives rise to the observed anti-correlation. We have placed 
HS2 assuming it was triggered by a disturbance induced by the huge initial 
flare peak at t w 28 minutes. The time-evolution of the H7 and BaC fluxes 
in Figure Id of KlOa indicates an apparent lack of new hydrogen Balmcr 
line emitting-regions during the secondary U-band peaks (we cannot defini- 
tively determine whether there is a newly-formed hydrogen Balmcr-cmitting 
(chromospheric ribbon) component cospatial with the hotspot because the 
observations are unresolved). Therefore, the disturbance likely propagated 
through the lower atmosphere, below the height of hydrogen Balmer line 
formation (upper chromosphere; J. Alfred, private communication 2010). 
Using a range of sound speeds in the lower atmosphere for the speed of the 
disturbance (« 5 — 10 km s^ 1 ), we find that HS2 is located at a distance 
that is approximately 30 — 60 Mm (Ryz CMi ~ ^00 Mm) from the site of 
the initial flare event Q 



The composite graphical model is preliminary (see Section [5]) and requires 
comparison to other complex flare events on dMe stars but especially to solar 
flares where we can spatially resolve each continuum component. Our group is 
currently working to obtain solar flare data that can be used to test the YZ 
CMi flare model using DST/ROSA and employing custom continuum filters 



(jJess et all l2010bl ; iKowalski et all 120111 ) 



4.1. Are the Secondary Flares Stellar Ellerman Bombs? 



Ellerman bombs are transient, compact pheno mena observed near evolving or 
emerging magnet i c field s in solar active regions ( Ellermanl . Il917 : Severnv . 1968 : 
Georgoulis et all . 120021 . and references therein). A typical signature of Eller- 
man bombs is emission in the wings and absorption in t h e core of Hct r e lative 
to the nea r by pl age intensity (e.g., Koval and Severnvl . Il970l : Bruzek . 1972 : 



Fang et all 120061) . The time-evolution prope rties include mean lifetimes of 



10 — 20 minutes and fine-structure variations ( Kurokawa et al. . 19821 : Qiu et al 



2000l ). In contrast to typical white-light flares, E llerman bomb s have symmetric 
light curves with similar rise and decay times ( Pavnel 1993 : Qiu et al. . 2000l : 



1 The white-light and hard X-ray footpoints have been observed to propagate along the 
polarity inversion lin e durin g large solar flares such a s the f amous 14 July 2000 flare 
^Fletcher and Hudson! J2001D ; iKosovichev and Zharkoval J2001I) ; lOiu et al\ J2010T ), but the 
spatial location of these kernels appears to change much faster, Ri 170 — 200 km s . 
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Figure 3. Graphic with continuum components (areas from t PS 126.5 minutes to scale) 
as they might appear in a spatially resolved observation. The BaC emitting region (yellow) 
resembles a two-ribbon structure (shown here as symmetric for simplicity) in a thin layer of the 
heated mid-to-upper chromosphere of previously reconnected magnetic loops. HS1 is shown 
as a collection of previously formed hot spots (purple) and HS2 is the proposed newly-formed 
continuum emitting region (white). The triangle indicates the assumed location of the inital 
flare peak, which generated a disturbance in the lower atmosphere that propagated into the 
surrounding active region and triggered the hot spots. The black circle helps orient the reader 
to the center of the star, which has a radius 0.3-Rq. The flare region is placed at an arbitrary 
location on the surfa ce. Several aspects of t his ca rtoon were inspired by observations of solar 
flare arcades, such as Fletcher and Hudson (2001). 



Jess et all l2010al ) . The Ellerman-bomb mechanism is not fully understood but 



has been attributed to magnetic reconnection in the low chromosphere (e.i 
Georeoulis eTaZl . l200zt . 



The secondary flares during the YZ CMi mcgaflarc exhibit several similari- 
ties to Ellerman-bomb phenomena on the Sun. Ellerman bombs have also been 
modeled phcnomcnologically as temperature bumps at or below the solar temper- 
ature minimum region (jFang et al. . 1200(1 iBerlicki. Heinzel. and Avretll2lnol) . 
The secondary YZ CMi flares have longer rise times (« two -five minutes) and 
are much more symmetric about the peak compared to other white light flares 
with similar total energy on YZ CMi (At r j sc « 0.5 — 1.8 minutes; MoffettJ . 



1974 : van den Oord et al . 1996). The absorption features of the blackbody- 



like continuum component are similar to the line-center absorption observed 
in Ha and Ca n during Ellerman bombs; unfortunately, our observations do not 
have sufficient spectral resolution to separate line-center and wing profiles. The 
preliminary finding (Section [4]) that the blackbody-like continuum component 
does not contain hydrogen-line emission may be consistent with magnetic re- 
connection taking place in the low atmosphere. However, in contrast to solar 



Ellerman bombs, which have been observed as a microflare trigger (|Jess et al 



2010a)), the secondary flares are possibly a consequence of the enormous YZ CMi 



flare peak event. 
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5. Summary and Future Work 



The time-resolved continuum data obtained during a megaflare on the dM4.5e 
star YZ CMi demonstrate the power of broad- wavelength coverage, low- resolution! 
spectra, which are unfortunately not currently available in the optical/near-UV 
for solar flares. In this manuscript, we show that the blackbody-like component 
(hot star-like emission) of the white-light continuum dominates the spectra dur- 
ing the secondary flares while the BaC (and Balmer lines) become less important; 
this observed anti-correlation is explained as continuum veiling. We combine the 
phenomenological m odels of the blackbody-likc component from KlOb and the 
Allred et al. ( 20061 ) RHD model spectrum of the BaC to reproduce the total 



flare emission at two times during the flare. The filling factors for the decaying 
BaC(+backwarming) component, a previously heated hot-spot component, and 
a newly heated hot spot component arc in the ratio of ~ 10 : 1 : 0.5. These areal 
coverages allow a comparison of each component to be made with observed solar- 
flare structures in large eruptive flare arcades. Although we generally assume that 
stellar-flare phenomena are simply "scaled- up" versions of solar-flare phenomena, 
one should not exclude the possibility that stellar flares might have fundamental 
differences in the white-light continuum, as the energies and timescales of dMe 
flares can be orders of magnitude larger than solar flares. New solar-flare obser- 
vations are needed to test the existence of the blackbody-like component and to 
better understand the properties of the BaC, which could be fully "unveiled" 
using spatially resolved solar observations. 

We provide evidence that the blackbody-likc component has several similar 

properties to solar El lerman bombs . A few solar flares have been known to exhibit 

the sp ectral features ( Svestka . 1963 ) and velocity characteristics ( Kosovichev and Zharkovai 



20011 ) of Ellerman bombs. Additional intensity-calibrated continuum measure- 



ments of Ellerman bombs and white-light flare kernels on the Sun, such as 
with DST/ROSA, would help illuminate the differences between these events 
and provide a comparison to stellar spectra of typical white-light flares and 
mcgaflarc-size events. 

Several aspects of the phenomenological models presented in this work are 
being improved. In addition to modeling hydrogen with more levels and includ- 
ing metallic transitions and molecular species, a more accurate consideration 
of charge balance is underway. The correct treatment of charge balance in a 
modified atmosphere is complicated by non-LTE ionization, but a new version 
of RH has been provided by H. Uitenbroek to account for this; the authors are 
currently working on a new suite of hot-spot models. 
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